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^ Abstract 

We consider the semiclassical rigid-body quantization of Skyrmion solutions of mass 
' numbers B — 4, 6, 8, 10 and 12. We determine the allowed quantum states for each 

Skyrmion, and find that they often match the observed states of nuclei. The spin and 
isospin inertia tensors of these Skyrmions are accurately calculated for the first time, and 
. are used to determine the excitation energies of the quantum states. We calculate the energy 

' level splittings, using a suitably chosen parameter set for each mass number. We find good 

(y-^ , qualitative and encouraging quantitative agreement with experiment. In particular, the 

' rotational bands of beryllium-8 and carbon-12, along with isospin 1 triplets and isospin 2 

On . quintets, are especially well reproduced. We also predict the existence of states that have 

' not yet been observed, and make predictions for the unknown quantum numbers of some 

observed states. 

^- 

■ 1 Introduction 

The SU(2) Skyrme model provides a geometrical picture of nuclear physics in which nuclei 
are identified with the topological soliton solutions of the model, known as Skyrmions [22] . 
The model has several advantages over conventional nuclear models. Firstly, single Skyrmions, 
which are identified with nucleons, are found to merge to some extent and to lose their indi- 
vidual identities in the solutions describing larger nuclei. This captures an important feature 
of nuclei with individual nucleons close together, and is something that can not be achieved 
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in conventional point nucleon models. Secondly, the Skyrme Lagrangian is defined in terms of 
only three parameters. For each mass number B we fit these parameters to the mass and size 
of the nucleus of zero isospin that has mass number B, to obtain predictions with reasonable 
quantitative accuracy. This is unlike potential and shell model calculations which require many 
finely tuned parameters. 

The model is a pion field theory, and is an extension of the nonlinear sigma model. The 
pion fields 7v{x) are combined into an SU(2)-valued scalar field, the Skyrme field 

U{x) = a{x)l2 + i7v{x) - r , (1) 

where r denotes the triplet of Pauli matrices and o"^ + tt • tt = 1. The Lagrangian density is 

m 

C = ^ Tr + ^ Tr [d^UU\d,UU^] [d^'UUl d'UU^ + ^mlF^ Tr {U - I2) , (2) 

where F.^ is the pion decay constant, e is a dimensionless parameter and is the pion mass. 
Using energy and length units of F^^/Ae and 2/e-F^ respectively, the Skyrme Lagrangian can be 
rewritten as 

L = y |-^TV(i?^i?'^) + ^Tr([i?,„i?,][i2^,i?n) + "^'Tr(?7-l2)|d=^x, (3) 

where i?^ = d^UW , and the dimensionless pion mass m = 2mT^/eFT^. As usual, the Lagrangian 
splits into kinetic and potential parts as L = T — V , with T quadratic in the time derivative of 
the Skyrme field. The potential energy V is given by 



V = / {-^ Tr (RiRi) - ^ Tr {[Ri, Rj][Ri, R^]) - m^TT{U - h)^ d'x . 



(4) 



The Skyrmions are minima of the potential energy and are static. They are labelled by their 
topological charge, B, the degree of the map C7 : — > SU(2), which is well-defined provided 
[/ — > I2 at spatial infinity, and is given by the integral 



B = j Bo{x)d-'x, (5) 

where 

B^^ix) = ^ e^^^p Tr [d'^UU^ d^UU^ df'UU^) . (6) 

We denote the minimised potential energy hy Mb- One interprets a charge B Skyrmion, after 
quantization, as a nucleus of mass number B. In this picture, nucleons and nuclei arise purely 
from the pion field and no explicit nucleonic sources are needed. 

It is necessary to semiclassically quantize the Skyrmions as rigid bodies. The Skyrme La- 
grangian is invariant under rotations in space. It is also invariant under the transformations 
U AUA^, where A G SU(2). This is isospin symmetry. The rotational and isorotational 
degrees of freedom are treated as collective coordinates and the Skyrmions are quantized in 
their rest frame by canonical methods. In this way the Skyrmions acquire spin and isospin. 
An advantage of this model over other nuclear models that involve collective rotational motion 
is that it incorporates isospin excitations. In the Skyrme model, the vacuum solution [/ = I2 
is isospin invariant, but for each classical Skyrmion solution, isospin symmetry as well as ro- 
tational symmetry is spontaneously broken. These symmetries are restored by the collective 
coordinate quantization. 
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The kinetic energy of a rigidly rotating Skyrmion (ignoring the rather trivial translational 
motion) is of the form 

T = -UiUijaj - ttiWijbj + -biVijbj , (7) 

where 5j and a,, are the angular velocities in space and isospace respectively, and Uij, Vij and 
Wij are inertia tensors |10l I16j . The momenta conjugate to bi and Oj are the body-fixed spin 
and isospin operators Li and Ki. The quantum Hamiltonian H is obtained by re-expressing T 
in terms of these operators. The space-fixed spin and isospin operators are denoted Jj and Ij 
respectively. Note that = and = K^. 

Finkelstein and Rubinstein showed that it is possible to quantize a B = 1 Skyrmion as a 
fermion, and showed that for even (odd) B the spin and isospin are integers (half-integers) 
m. Discrete symmetries of the classical Skyrmion solution^ give rise to further Finkelstein- 
Rubinstein (FR) constraints on the space of quantum states |^'). These constraints are of the 
form 

gi92n2-Lgieini.K|^^ = XFr|^> , (8) 

where ni, n2 and 6i, 62 are, respectively, the axes and angles defining the rotations in isospace 
and space associated with a particular symmetry, and xfr = ±1- The FR signs, xfr, define 
a one-dimensional representation of the symmetry group of the Skyrmion. Krusch has found 
a convenient way to calculate them [13], and we use this method here. A basis for the wave 
functions |^') is given by the products | J, L3) ® |/, ^^3), the tensor products of states for a rigid 
body in space and a rigid body in isospace. J and I are the total spin and isospin quantum 
numbers, L3 and the projections on to the third body-fixed axes, and the space projection 
labels (which are the physical third components of spin and isospin) are suppressed. The 
FR constraints only allow a subset of these products as physical states. A parity operator is 
introduced by considering a Skyrmion's reflection symmetries. Quantum states are therefore 
labelled by the usual quantum numbers: spin-parity J'^ and isospin /. 

The inclusion of the third term in the Lagrangian density, which involves the pion mass, 
has a significant effect on the shapes and symmetries of the classical Skyrmion solutions. This 
effect is more marked for larger values of B. For zero pion mass, the Skyrmions with B up to 
22 and beyond resemble hollow polyhedra, with their baryon density concentrated in a shell of 
roughly constant thickness, surrounding a region in which the baryon density is very small [8]. 
This disagrees with the approximately uniform baryon density observed in the interior of real 
nuclei. Fortunately, it has been established that the hollow polyhedral solutions for i? > 8 do 
not remain stable when the pion mass is set at a physically reasonable value, with m ~ 1 [9]. 
This is because in the interior of the hollow polyhedra the Skyrme field is very close toU = — 12, 
and here the pion mass term gives the field a maximal potential energy, and hence instability. 
This instability results in the interior region splitting into separate smaller subregions. The 
stable Skyrmion solutions are found to exhibit clustering: small Skyrmion solutions, such as 
the cubically symmetric B = 4 solution, appear as substructures within larger solutions [7j. This 
is a very encouraging development as it has been believed for some time that alpha-particles 
exist as stable substructures inside heavier nuclei. The most remarkable success of the alpha- 
particle model is in its prediction for the binding energies of nuclei which can be formed out of 
an integral number of alpha-particles. 

In earlier work, the B = 4, 6 and 8 Skyrmions have been quantized [19^ I16j . using some 
approximations to the Skyrmion solutions based on the rational map ansatz [12]. The inertia 
tensors used had the right symmetries, but not the correct numerical values. In this paper 

^Only the B = 1 and B = 2 Skyrmions possess continuous symmetries. 
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we consider the S = 4, 6, 8, 10 and 12 Skyrmions, using a consistent numerical scheme to 
recalculate all these solutions. Each of these Skyrmions can be viewed as being built up from 
B = 4 cubes (it is possible to regard the B = 6 Skyrmion as being made up of a -B = 4 cube 
and a, B = 2 torus, and the i? = 10 solution as consisting of two B = A cubes together with 
two B = 1 Skyrmions). We numerically relax field configurations to the stable Skyrmions, for 
various values of the dimensionless pion mass m, and then compute their static energies Mb, 

1 /2 

charge radii {^r^) , and inertia tensors. Appendix A tabulates the calculated numerical values 
for m = 0.5, 1 and 1.5, and describes the quadratic interpolation method used to estimate these 
quantities for any given m. In the next section we describe our method of choosing m and of 
calibrating the model to ensure that it provides quantitatively accurate predictions of nuclear 
properties. 

For each Skyrmion we determine all the FR-allowed quantum states and their excitation 
energies, working up to spin and isospin values just beyond the edge of what is experimentally 
accessible. Our rigid-body quantization leads to an infinite tower of quantum states (the B = 1 
Skyrmion, for example, has quantum states for all half-integer values of spin). However, in 
practice we expect a Skyrmion to deform as it spins (this is known to occur for the B = 1 
Skyrmion [6]), and this disallows many higher spin states. We are therefore not concerned that 
in some cases we predict higher spin states that are not experimentally observed. It is also 
possible that a Skyrmion might break up as it spins, a phenomenon that is known to occur for 
real nuclei. 

2 Calibration 

The model was first calibrated by Adkins, Nappi and Witten, by fitting the model in the B = 1 
sector to the masses of the proton and delta resonance [21 tt] • In a recent paper a new calibration 
was obtained by equating the mass and size of the quantized B = 6 Skyrmion to the mass and 
size of the lithium-6 nucleus [19] . However this was performed using the approximate Skyrmion 
found using the rational map ansatz. Evidently, there are many possible ways in which the 
model may be calibrated. 

The three parameters of the model are the pion decay constant (experimentally 186 MeV), 
the pion mass rriT^ (experimentally 138 MeV), and the dimensionless parameter e. Strictly speak- 
ing, one may argue that we are only free to set e, as the other constants are fixed by experiment. 
However, we consider it permissible to vary F^r, as we consider it to be a renormalised pion de- 
cay constant. We also allow e to vary with mass number, e = e{B), in order to get better 
agreement with experiment. The length and classical energy scales are 21 eF^^ and F-j^/Ae^ re- 
spectively. We recall that the quantum Hamiltonian for a rigidly rotating body is equal to 
the squared angular momentum operator divided by twice the moment of inertia of the body 
|15] . The moment of inertia has units equal to the mass scale multiplied by the square of the 
length scale: {F-j^/Ae) x {2/eFT^Y = l/e^-^Tr- The quantum energy scale is its reciprocal, e^F^^. 
The total energy of a quantum state of a Skyrmion is therefore equal to {F-,^/4:e)M. b + ^^F-kEq, 
where Eq is the quantum kinetic energy of the state. The ratio of the total energy of a quantum 
state to the classical energy of the Skyrmion is 1 + Ae^{EQ/ Mb)- The ratios of the quantum 
energies of Skyrmion states are therefore sensitive to the value of e. 

The classical Skyrmion solutions match the experimental pion tails of nuclei if we use the 
physical value of m-Tr of 138 MeV. For this reason, we keep rnvr fixed at this value. For each B, we 
choose m (and therefore fix the length scale, as m = 2m-,^/eF-,^) such that the calculated mean 
charge radius agrees with that of the nucleus with zero isospin with this value of B. Within 
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the Skyrme model, the mean charge radius of a nucleus with zero isospin is estimated to be the 
square root of 

2^ / Bq{x) (fx 



f Bo{x) (Pi 



(9) 



since the electric charge density is half the baryon density |19j . We then set the classical energy 
scale such that the Skyrmion mass agrees with the nuclear mass, which is done by setting 
F-j^/Ae to be the nuclear mass divided hy Mb- We use the actual nuclear masses, not taking 
into account small quantum energies which are of the order of 0.1% of the nuclear mas^. It is 
found that in each case m takes a value between 0.6 and 1.2. The experimental data, the values 
of m used, and the length and energy scales are listed in Table lH. 
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Nucleus 


Mean charge 


Mass 


m 


Length 


Classical 


Quantum 






radius 


(MeV) 




scale 


energy scale 


energy scale 






(fm) 






(fm) 


F^/Ae (MeV) 


e^F^ (MeV) 


4 


^He 


1.71 


3727 


0.820 


1.173 


6.168 


4589 


6 


6Li 


2.55 


5601 


1.153 


1.648 


5.753 


2491 


8 


8Be 


2.51 


7455 


0.832 


1.190 


6.336 


4339 


10 


lOg 


2.58 


9327 


0.830 


1.187 


6.354 


4350 


12 


12c 


2.46 


11178 


0.685 


0.980 


6.527 


6214 



Table 1: Experimental data and calibration for each B. 



A larger length scale, for example in the case of = 6, takes into account loose vibrational 
motion, and leads to larger moments of inertia, which is desirable. The small length scale for 
i? = 12 takes into account the compact size of the carbon-12 nucleus. Another reason to use 
separate parameter sets for each B comes from comparing the rotational bands of beryllium-8 
and carbon-12. These nuclei have O"*" ground states and 2"*" and 4+ excited states, at 3.0 MeV and 
11.4 MeV respectively for beryllium-8, and at 4.4 MeV and 14.4 MeV respectively for carbon-12. 
Naively, one might expect carbon-12 to be larger and heavier than beryllium-8 (in dimensionless 
units, the B = 12 Skyrmion is larger in size and has a larger classical mass than the 5 = 8 
Skyrmion), and for it therefore to have a larger moment of inertia. As the moment of inertia 
appears in the denominator of the quantum Hamiltonian, this would lead to carbon-12 having 
a smaller rotational band splitting than that of beryllium-8, which is not the case. Evidently, 
the only way to deal with this problem in our model is to use a different parameter set for 
each. Experimentally, the mean charge radii of the nuclei we are considering are approximately 
constant for Q < B < 12, whereas the dimensionless Skyrmion mean charge radius increases with 
B. The optimised parameter sets allow us to keep the nuclear mean charge radii approximately 
constant. 

^ Lithium-6 and boron-10 have quantum energies associated with the non-zero spins of their ground states (1^ 
and S"*" respectively). We calculate the energy of the l""" state of the 5 = 6 Skyrmion to be 1.6 MeV (see section 
4), which is negligible compared to the mass of lithium-6, 5601 MeV. An estimate of the boron-10 mass with its 
quantum energy removed can be obtained by taking the average of the masses of beryllium-10 and carbon-10 
and subtracting from this our estimate of the O""" state with isospin 1, 10.4 MeV (see section 6). This yields a 
value of 9319 MeV, which is very close to the accepted mass of boron-10, 9327 MeV. 

■^The mean charge radius of beryllium-8 has not been measured, due to its instability. Here we use the value 
for its isobar lithium-8. 
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3 5 = 4 



The minimal-energy B = 4 Skyrmion has octahedral symmetry and a cubic shape. A surface 
of constant baryon density is presented in Fig. 1. The colour scheme represents the direction 
in isospace of the associated pion fields. For regions in space where at least one of the pion 
fields does not vanish, the normalised pion field tt can be defined, and takes values in the unit 
sphere. We colour this sphere by making a region close to the north pole white and a region 
close to the south pole black. On an equatorial band, where tts is small, we divide the sphere 
into three segments and colour these red, blue and green. In Fig. 1, opposite faces share the 
same colour and vertices alternate between black and white. 

For a derivation of the Skyrmion's quantum states we refer the reader to Ref. p6], and 
here we only state the main results. For an earlier discussion, see Ref. |28j . The rotational 
symmetry group of the Skyrmion, Oh, is one of the largest point symmetry groups. This leads 
to particularly restrictive FR constraints on the space of allowed states. The Oh symmetry 
implies that the inertia tensors are diagonal with Un = U22 and C/33 different, Vij proportional 
to the identity matrix and Wij = 0. The quantum collective coordinate Hamiltonian is therefore 
the sum of a spherical top in space and a symmetric top in isospace, 

H = + J—I^ +(— —] Kl . (10) 

The inertia tensors are given in Appendix A.l. The lowest state is a 0^ state with isospin 0, 
which has the quantum numbers of the alpha-particle in its ground state. The first excited 
state with isospin is a 4"*" state, which has not been experimentally observed, because of 
its high energy. The lowest state with isospin 1 is a 2~ state, which matches the observed 
isotriplet including the hydrogen-4 and lithium-4 ground states. The energy levels are given 
in Table 2. Only the quantized kinetic contributions to the total energies are listed; the static 
Skyrmion mass must be added to give the total energy. For example, the ground state has zero 
kinetic energy and therefore its total energy is precisely 3727 MeV, as given in Table 1. The 
4"*" state has additional kinetic energy of IO/V33 = 129.7 x 10~^ in dimensionless units. In the 
final column of Table 2 we list the values in physical units, using the i? = 4 conversion factor 
c^Ft^ = 4589 MeV given in Table 1. We overpredict the excitation energy of the 2~ isotriplet as 
43.2 MeV compared with an average experimental value of 23.7 MeV [21] . 

In summary: The ground state of helium-4 and the isotriplet of 2~ states arise as quantum 
states of the = 4 Skyrmion. The cubic shape of the Skyrmion may become deformed as it 
spins. It may be for this reason that the 4"^ state is not observed experimentally. 




Figure 1: A surface of constant baryon density for the B 
indicate different directions of the pion fields. 



4 Skyrmion. Different colours 
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i^ivie V ) 
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4+ 


129.7 


59.5 
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2~ 


94.1 


43.2 



Table 2: Energy levels of the quantized B = 4 Skyrmion. 

4 5 = 6 

The B = 6 Skyrmion has symmetry (see Fig. 2). We refer the reader to Ref. [16j for 
a discussion of its quantization. The quantum Hamiltonian is that of a system of coupled 
symmetric tops: 

h = -^j^ + -^i^+(-^-^)lI + (^-^)kI + ^l,Ks, (11) 

2Vn ^ 2Uu V2A33 2VnJ ' V2A33 2?7ii ^ ^ ^ A33 ' ' ' ^ ^ 

where A33 = C/33V33 — W^^. Its allowed quantum states are listed in Table 3, together with 
their energy levels computed using the inertia tensors in Appendix A. 2. 




Figure 2: A surface of constant baryon density for the B = Q Skyrmion (two viewpoints). 



The Skyrme model qualitatively reproduces the experimental spectrum of lithium-6 and 
its isobars, and predicts some further states, including J'^ = 4~, 5"*" and 5~ excited states 
of lithium-6 with isospin 0. The ground state of lithium-6 is correctly predicted to be a 1^ 
state. We also find a 3^ excited state, which is seen experimentally. However we overpredict 
its excitation energy by roughly a factor of five. The model does not account for centrifugal 
stretching nor the allowed decay of the 3"^ state to an alpha-particle plus a deuteron, which 
may be the reason for our overprediction. The lowest allowed state with isospin 1 is a 0^ 
state, which is seen experimentally as an isotriplet which includes the helium-6 and beryllium- 
6 ground states. An excited 2+ state of this isotriplet exists, which we also find within our 
model. However again we overpredict its excitation energy. We predict an additional 2^ state 
with isospin 1. The lowest allowed negative parity state with isospin 1 is predicted to be a 2~ 
state with excitation energy 18.2 MeV. We therefore suggest that the observed 9.7 MeV state 
of helium-6 is our second 2+ state. Lithium-6 has a 2~ state with isospin 1 at 18.0 MeV. We 
predict that the 18.7 MeV state of helium-6 has J'^ = 2~ , and is one of its isotriplet partners. A 
2~ state of beryllium-6 is observed at 29.1 MeV. Perhaps this state completes the isotriplet, but 
its high energy makes this unclear. States of lithium-6 and beryllium-6 with J'^ = 3~ and 4~ 
and with isospin 1 have been experimentally observed. We predict these states, with roughly 
the correct excitation energies. We also predict the existence of 3^ and 4"'" states with isospin 
1, which have not been seen experimentally. The ground state of hydrogen-6 has isospin 2, 
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on 1 
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34.2 






1 -i r\\ .- — 1 1 c\ r\\ 

|1,0) (g |2,0) 
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37.0 






|1,0) ® (|2,2) + |2, -2)) 
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35.8 




9 + 


lo o\ |0 1\ |0 0\ fO^ |0 1\ 


10 / .o 








|2,0)®(|2,2) -|2,-2)) 


156.9 


39.1 






|2,2)«)|2,-1)-|2,-2)«)|2, 1) 


167.8 


41.8 



Table 3: Energy levels of the quantized S = 6 Skyrmion. 



at 28.2 MeV above the lithium-6 ground state, and an undetermined spin. The Skyrme model 
gives the lowest allowed state with isospin 2 as a 0~ state, and its excitation energy to be 
34.2 MeV. Higher spin excited isospin 2 states are also allowed in the model, but they have not 
been experimentally observed. 

In summary: The quantum numbers of the low-lying states of the B = 6 Skyrmion agree 
with those of the nuclei of mass number 6, although we overpredict their excitation energies. Our 
calculated values for the excitation energies of the isospin 1 states are, however, quantitatively 
good. We have also made predictions for the spins of two excited states of helium-6, and predict 
that the hydrogen-6 ground state is a 0"" state. 
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4.1MeV J=0+I=1 

Helium-6 
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3.1McV J=OtI=l 

Beryllium-6 



Figure 3: Energy level diagram for nuclei of mass number 6. Mass splittings between nuclei 
are adjusted to eliminate the proton/neutron mass difference and remove Coulomb effects, as 
described in Ref. l25J. 



5 B = 8 



When m = 1, the stable B = 8 Skyrmion is D4/i-symmetric, and resembles two touching 
B = 4 cubes, matching the alpha-particle model picture of beryllium-8 as an almost bound 
configuration of two alpha-particles. A surface of constant baryon density is displayed in Fig. 
4. In Ref. we approximated the B = 8 Skyrmion as a "double cube" of two B = A 

Skyrmions, and made estimates for its inertia tensors, enabling us to estimate its energy levels. 
We refer the reader to this paper for a discussion of its quantization. The quantum Hamiltonian 
is the sum of a symmetric top in space and an asymmetric top in isospace: 
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2Vu 



J2 + 



1 
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2^33 2^11 



1 



2Ur 



-Kl 



1 



2C/22 



Kl + 



1 



-Ki 



(12) 



The exact inertia tensors are given in Appendix A. 3. As anticipated on symmetry grounds, the 
inertia tensor Uij has three distinct eigenvalues (the earlier double cube estimate made two of 
them equal). The energy levels calculated using the exact inertia tensors are listed in Table 4. 

Figure 4 is an energy level diagram for nuclei of mass number 8. The Skyrme model predic- 
tions for positive parity states agree well with experiment. The ground state of beryllium-8 is 
correctly determined to be a O"*" state. The rotational band of beryllium-8 is remarkably well 
reproduced in our model: we predict the 2"*" and 4"*" states at 3.0 MeV and 10.2 MeV respec- 
tively, which is very close to the experimental values of 3.0 MeV and 11.4 MeV respectively. 
We predict a second 4"*" state with isospin at 24.3 MeV. Experimentally, beryllium-8 has two 
further 4+ states with isospin 0, at 19.9 MeV and 25.5 MeV. 

The model predicts that the lowest allowed state with isospin 1 has = 0~. However, the 
ground states of lithium-8 and boron-8 are believed to have J'^ = 2^. Perhaps the 0~ isotriplet 
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Figure 4: A surface of constant baryon density for the B = 8 Skyrmion, resembling two B = 4 
cubes. 
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23.2, 27.2, 37.4 


2 


0+ 


2 


87.6, 93.5 


38.0, 40.6 




0" 


1 


90.9 


39.4 




2+ 


3 


94.6, 100.5, 108.1 


41.0, 43.6, 46.9 




2- 


2 


98.0, 103.0 


42.5, 44.7 



Table 4: Energy levels of the quantized B = 8 Skyrmion. n is the number of FR-allowed states 
with given / and J'^ . 

may be seen experimentally in the future, for it is known that low-lying spin 0, negative parity 
states are difficult to observe, as experienced in the search for the bottomonium and charmonium 
ground state mesons r]ii and rjc [5, 21j. Our estimate for the excitation energy of the isotriplet 
is 19.3 MeV, which experimentally has an average excitation energy of 16.7 MeV. In addition, 
we predict a 3^ isotriplet at 22.4 MeV, to be compared to the experimental value of 19.0 MeV. 
The model forbids spin 1 states with isospin 1. However several 1"*" states with isospin 1 have 
been observed in the spectrum of lithium-8. Perhaps these may arise as quantum states of 
an alternative B = 8 Skyrmion, or from the quantization of further degrees of freedom, such 
as vibrational modes. We predict two 2~ states with isospin 1, at 16.1 MeV and 20.1 MeV. A 
state with these quantum numbers has been seen in the spectrum of boron-8, at 19.9 MeV. We 
predict a 4"'" state with isospin 1 at 26.5 MeV, which is seen in the spectrum of lithium-8 at 
23.5 MeV, although its beryllium-8 and boron-8 partners are not yet confirmed. Beryllium-8 
has a 4~ state at 20.9 MeV, which we identify with our 4~ state with isospin 1 at 23.2 MeV. 

We find that the lowest state with isospin 2 has J'^ = 0"*". Experimentally this forms an 
isospin 2 quintet which includes the ground states of helium-8 and carbon-8. We calculate its 
excitation energy to be 38.0 MeV, to be compared to the experimental value of 27.3 MeV. A 
state of helium-8 has been observed at approximately 31.1 MeV, which we predict at 41.0 MeV. 

In summary: The spectrum of beryllium-8 is very well reproduced in our model, in particular 
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-31.1 MeV J=2+ 



28.1 MeV J=0+. 1=2 



Helium-8 



22.4 MeV J=l+, 1=1 - 



17.0 MeV J=l+, 1=1 - 



27.8 MeV J=0+, 1=2 



23.5 MeV J=4+. 1=1 



20.2 MeV J=l+, 1=1 



19.3 MeV J=3+, 1=1 



27.5 MeV J=0+, 1=2 



25.5 MeV J=4+, 1=0 



24.0 MeV 



J=(l,2)7l=0 
20.9 MeV 1=4- 



19 9 MeV 1=4*. 1=0 
19.1 MeV 1=3+, 1=1 



16.6 MeV 1=2+, 1=1 



17.0 MeV J=2+. 1=1 

Lithium-8 



11.4 MeV J=4+, 1=0 



27.0 MeV 1=0+, 1=2 



19.9 MeV 1=2-, 1=1 



18.6 MeV 1=3+, 1=1 



26.3 MeV J=0+, 1=2 



Carbon-8 



17.1 MeV J=l+, 1=1 



16.4 MeV 1=2+, 1=1 

Boron-8 



3.0 MeV I=2+, 1=0 



J=0+, 1=0 

Beryllium- 8 

Figure 5: Energy level diagram for nuclei of mass number 8. Mass splittings between nuclei 
are adjusted to eliminate the proton/neutron mass difference and remove Coulomb effects, as 
described in Ref. 



its isospin rotational band. So too are the experimentally observed isospin 1 triplets and the 
isospin 2 quintet. Predicted 0~ states with isospin 1 have not been seen in the spectra of nuclei 
of mass number 8, however as we mentioned they may be difficult to observe. We have been 
unable to explain the appearance of 1^ states of lithium-8 and boron-8. We are led to consider 
refinements of our model and its quantization in order to address this. 



6 B = 10 

In Ref. [9\ the minimal-energy solution for B = 10 was found to have D2h symmetry. This 
Skyrmion can be thought of within the context of the alpha-particle model as a pair of B = 4 
cubes with two B = 1 Skyrmions between, as illustrated in Fig. 6. 

Here, we quantize this Skyrmion for the first time. We use the rational map ansatz to de- 
termine its FR constraints. For an overview of the ansatz, which provides close approximations 
to the exact Skyrmion solutions, see Ref. p^. While it does not provide quantitatively exact 
results, it precisely describes the symmetry group of many Skyrmions, and therefore can be 
used to derive their FR constraints. In the B = 10 case a suitable rational map is 

, . _ g + + + 6z^ + ez^ + 
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Figure 6: A surface of constant baryon density for the B = 10 Skyrmion (two viewpoints) 



with a = 0.28, /? = -9.37, 7 = 14.83, 5 = 4.98 and e = 3.02. The D2 rotation group is 
generated by two C2 symmetries, which correspond to the following symmetries of the rational 
map: 

R{-z) = R{z), R{l/z) = l/R{z). (14) 
The associated FR constraints are determined using Krusch's method |14] to be 

e^'^^-^\^) = \^) , e''^^^e''^^^\^) = -\^) , (15) 

so the signs xfr form one of the non-trivial one-dimensional representations of D2 ■ 

Generally, the parity operator in the Skyrme model is an inversion in space and isospace: 
V : C/(x) ?7^(-x). A rational map is invariant under parity if it satisfies 



R{-l/z) = -1/R{z) . (16) 

For a quantized Skyrmion described by such a rational map, all states would have positive 
parity. However, as the rational map (I13p satisfies the inversion symmetry 



R{-l/z) = l/Riz) , (17) 

the parity operation in this case is equivalent to a single rotation in isospace, given hyV = e*'^^^ . 
The parity of a quantum state is its eigenvalue when acted upon by this V. Note that we attach 
the parity label to the spin quantum number, to form J'^ , as is conventionally done in nuclear 
physics, despite the fact that in the Skyrme model the parity operator generally reduces to a 
combination of rotations in space and in isospace. 

For the B = 10 Skyrmion, the symmetries imply that the inertia tensors Uij and Vij are 
diagonal, and the only non-zero component of the mixed inertia tensor Wij is W33 . The quantum 
Hamiltonian is that of a system of coupled asymmetric tops: 

2^11 ' 2V22 ' 2A33 ' 2C/n ' 2U22 ' 2A33 3 ^ A33 ' ' ' ^ ' 

where A33 = J733 V33 — W33 as before. In the absence of symmetry and FR constraints, (2J -|- 
1) X {21 + 1) different non-degenerate levels would correspond to any given pair of J and I. 
Imposing the FR constraints, however, substantially reduces the number of allowed states. The 
calculation of energy levels is similar to the case of a general asymmetric top, described in 
Ref. |jl5j. However, the final term in ()18p mixes states of the form |J, L3) -|- |J, — -L3) and 
|J, L3) — |J, — L3) (and similarly for isospin basis states). In Table 5 we list the total number 
of allowed states, n, for different combinations of spin and isospin, together with their energy 
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eigenvalues. The energy levels are calculated by diagonalizing the Hamiltonian in matrix form, 
separately for each combination of J'^ and /. The precise forms of the eigenvectors are omitted 
as they do not add anything to our discussion. 

The ground state of boron-10 has J'^ = S"*" and its first excited state has J'^ = 1+ at 
0.7 MeV. We incorrectly determine the l"*" state to be the ground state, and 3^ states to be 
excited states. However, this problem arises in other models of boron-10, for example in models 
involving nucleon-nucleon potentials in chiral perturbation theory [20j. Boron-10 has further 
isospin excited states, including 2+, 3^ and 4^" states at 3.6 MeV, 4.8 MeV and 6.0 MeV 
respectively. We find that our model only allows positive parity states with isospin 0, and our 
predictions for the excitation energies of the aforementioned states are of the correct order of 
magnitude. Second 3"*" and 4"'" states with isospin are allowed, which we identify with the 
states of boron-10 at 7.0 MeV and 10.8 MeV, respectively. Curiously, 2~, 3~ and 4~ states of 
boron-10 with isospin have been found experimentally, at 5.1 MeV, 6.1 MeV and 6.6 MeV, 
respectively. Consideration of further degrees of freedom or the quantization of an alternative 
B = IQ Skyrmion may be necessary to account for these negative parity states. 

We find that the lowest allowed state with isospin 1 has J'^ = 0"^, and is at 4.9 MeV above 
the lowest 3^ isospin state. Experimentally this is observed as an isospin 1 triplet that 
includes the beryllium-10 and carbon-10 ground states, and has an average excitation energy 
of 1.8 MeV relative to the boron-10 ground state. Three 2+ spin excitations of this isotriplet 
have been observed, at average excitation energies 5.2 MeV, 7.4 MeV and 8.9 MeV. We predict 
two such excitations at 7.6 MeV and 8.3 MeV. Note that our model disallows a J'^ = l'^ state 
with isospin 1. This agrees with experimental observations. Additionally, 1~, 2~ and 3~ states 
with isospin 1 have been seen in the spectrum of beryllium-10, and our predictions for their 
excitation energies are close to the experimental values. A number of 2~ states with isospin 1 
are also seen in the spectrum of boron-10. 

For isospins 2 and 3, we find that all possible spins and parities are allowed, and that in 
both cases the lowest state is a 0+ state. The spins of the lithium-10 and nitrogen-10 states 
(both with isospin 2) are not well established. The lowest state of lithium-10 is at 23.3 MeV 
above that of boron-10, and its spin is uncertain: either 1~ or 2~. We calculate these states to 
exist at 26.7 MeV and 28.8 MeV above the lowest 3^ isospin state. Lithium-10 has an excited 
1^ state, and the lowest state of nitrogen-10 is believed to be a l"*" state. We predict two 1^ 
states with isospin 2, and slightly overpredict their excitation energies. The lowest helium-10 
state is a 0"*" state at 39.4 MeV above the boron-10 ground state. This is to be compared to 
our value of 62.2 MeV. In our model, the isobar splittings increase in proportion to /(/ + 1). 
However, as can be seen from the nuclear energy level diagrams, this approximately quadratic 
behaviour is not precisely reflected in the data. 

In summary: The quantum numbers and excitation energies of the states of nuclei of mass 
number 10 are reasonably well described by our model. However, the appearance of a 1^ 
state as the isospin ground state disagrees with that of boron-10, but as we mentioned this 
is a well-known artifact of nuclear models. The negative parity states with isospin in the 
nuclear spectra may arise as quantum states of an alternative i? = 10 Skyrmion, invariant 
under an alternative symmetry group which may permit such states. We find significantly more 
allowed states with isospins 2 and 3 than have been found experimentally. However we are not 
too concerned about those of higher spin, as the D2h symmetry may become deformed as the 
Skyrmion spins. 
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Figure 7: Energy level diagram for nuclei of mass number 10. Individual isobars are shifted 
vertically for clarity, and mass splittings between nuclei are adjusted to eliminate the pro- 
ton/neutron mass difference and remove Coulomb effects, as described in Ref. [26]. 
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/ 


J"" 


n 


E (xlO-*) 


E (MeV) 





1+ 


1 


2.5 


1.1 




2+ 


1 


6.6 


2.9 




3+ 


2 


12.7, 16.1 


5.5, 7.0 




4+ 


2 


21.1, 24.3 


9.2, 10.6 


1 


0+ 


1 


23.8 


10.4 




0" 


1 


24.9 


10.8 




1- 


1 


27.5 


11.9 




2+ 


2 


30.0, 31.7 


13.1, 13.8 




2" 


3 


31.1, 31.6, 32.8 


13.5, 13.8, 14.3 




3+ 


1 


37.8 


16.5 




3~ 


3 


37.6, 38.9, 41.1 


16.4, 16.9, 17.9 




4+ 


3 


44.4, 46.1, 51.1 


19.3, 20.1, 22.2 




4~ 


5 


45.5, 46.2, 47.2, 49.3, 


19.8, 20.1, 20.5, 21.5, 








52.2 


22.7 


2 


0+ 


1 


76.0 


33.1 




0~ 


1 


72.7 


31.6 




1+ 


2 


73.9, 78.5 


32.2, 34.2 




1" 


1 


74.0 


32.2 




2+ 


4 


78.1, 82.2, 82.7, 83.9 


34.0, 35.8, 36.0, 36.5 




2" 


3 


78.9, 79.2, 80.6 


34.3, 34.4, 35.0 




3+ 


5 


84.2, 87.6, 88.7, 90.1, 


36.6, 38.1, 38.6, 39.2, 








92.1 


40.1 




3" 


3 


85.2, 86.7, 88.6 


37.1, 37.7, 38.5 


3 


0+ 


2 


142.9, 147.5 


62.2, 64.2 






2 


144.1, 153.2 


62 7 66 7 

\J • 1 ^ \J\J m 1 




1+ 


1 


150.0 


65.2 




1" 


2 


146.8, 155.7 


63.9, 67.7 




2+ 


5 


149.1, 150.8, 153.7, 154.2, 


64.9, 65.6, 66.9, 67.1, 








155.3 


67.6 




2- 


6 


150.0, 151.0, 151.7, 159.3, 


65.3, 65.7, 66.0, 69.3, 








159.9, 161.1 


69.6, 70.1 




3+ 


4 


156.9, 160.2, 161.6, 163.6 


68.3, 69.7, 70.3, 71.2 




3- 


6 


157.0, 157.9, 160.4, 165.9, 


68.3, 68.7, 69.8, 72.2, 








167.2, 169.3 


72.8, 73.7 



Table 5: Energy levels of the quantized 5 = 10 Skyrmion. 
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7 B = 12 



In the alpha-particle model, the classical minimum of the potential energy for three alpha- 
particles occurs when they are arranged in an equilateral triangle. The minimal-energy solution 
of the Skyrme model in the B = 12 sector has Cs^ symmetry, but there is a solution of 
very slightly higher energy with a larger D^h symmetry, which is a saddle point and not a 
local energy minimum. Both have an equilateral triangle shape. Here, we quantize the D^/^- 
symmetric solution (see Fig. 8) , which we believe provides a physically more realistic picture of 
the nucleus. A more refined analysis might include an anharmonic vibrational mode centred on 
the L'a/i-symmetric solution and oscillating through two Ca^-symmetric minima. This solution 
can be approximated using the double rational map ansatz [17] . We use this to determine its 
FR constraints. The ansatz involves a Ds/j-symmetric outer map of degree 11, and a 

spherically symmetric degree 1 inner map, i?™, together with an overall radial profile function. 
The maps are [7]: 

^' z^-fz^+pz^ + az^ + l)' ^ ' 

Ri'^(z) = -z, (20) 

where a = —2.47, P = —0.84 and 7 = —0.13. The orientation of the inner map is chosen so 
that the D3/J symmetry is realised in a way compatible with the outer map. Both maps satisfy 

R{e'^ z) = e'^ R{z) , R{1/ z) = 1/ R{z) . (21) 

As the baryon number is a multiple of four, the FR signs form the trivial representation of -D3, 
and so the FR constraints are IT] 



gi^i3gi^i^3|^) = 1^) , e^^^i 6*^-^1 1*) = I*) . (22) 
Both maps satisfy the reflection symmetry 



Ril/z) = l/Riz) , (23) 

and so the parity operator is equivalent to V = e*'^^3g«'r^3 _ ^pj^g symmetry implies that 
the inertia tensors are diagonal, with Uu = U22, Vu = V22 and Wu = W22, so the quantum 
Hamiltonian is that of a system of coupled symmetric tops: 

where M = L + K, A33 is as before, and An = UuVu — W^^. 

The states that are allowed by the FR constraints are listed in Table 6. Each of the al- 
lowed isospin states is also an eigenstate of the Hamiltonian, and their quantum energies 
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Figure 8: A surface of constant baryon density for the Ds/j-symmetric B = 12 Skyrmion (two 
viewpoints). 



Ej^,i,\L3\,\K3\ are: 



-^0+, 0,0,0 
-^2+, 0,0,0 



0,3,0 
4- ,0,3,0 
-^4+ ,0,0,0 
-^5- ,0,3,0 
-^6^,0,3,0 
-^6+ ,0,0,0 
-^6+ ,0,6,0 



0, 

3f/ii/An , 

3C/ii/2An + 9C/33/2A33, 
llC/ii/2An + 9C/33/2A33: 
lOf/ii/An , 

21C/n/2An + 9?733/2A33: 
33C/ii/2An + 9C/33/2A33: 
21C/n/An , 

3C/ii/An + 18t/33/A33 • 



(25) 



These states also result from the rigid-body quantization of an equilateral triangle with alpha- 
particles at its vertices, and are not a prediction characteristic of the Skyrme model itself. The 
states fall into rotational bands labelled by IL3I = 0, 3, 6, .... As in Refs. [21 H], we suggest 
that the second experimental 2~ state of carbon-12 at 13.4 MeV has been misidentified, and is 
really a 4~ state. Again as in Ref. we predict a relatively low-energy 6"*" state of carbon-12, 
with IL3I =6. Such a state has not yet been seen experimentally. 

The isospin 1 states in Table 6 are not all individually eigenstates of the Hamiltonian, as 
they are not eigenstates of M'^. It is convenient to introduce basis states |J, /;M, M3), and 
express the states we have in terms of these. For example, for the two orthonormal FR-allowed 
1^ states: 



1+,1,0,0/ 



|1,0) |1,0) 



'3|i,i;2,o) 



^|i,i;0,o), 



(26) 



1+.1,1 



,1> 



V2 



ll,-l) + ll,-l) 



ll,l)) 



1,1;2,0) + W-|1,1;0,0).(27) 



Neither is an eigenstate of the Hamiltonian, as neither is an eigenstate of M^, but we can find 



two linear combinations of 



+ ,1,0,0/ 



and V^i 



+ ,1,1,1/ 



that are eigenstates by diagonalizing a 



matrix Hamiltonian. Their energies are the eigenvalues of the 2x2 matrix 



'^i+, 1,0,0 
.*i+,i,i,i 



H 
H 



*i+,i,o,o/ 
^i+,i,o,o/ 



(^1+, 1,0,0 

(^1+,1,1,1 



H 
H 



,1,1,1/ 
,1,1,1/ 



(28) 
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where 



a = C/ii/An + yn/An, (29) 

b = C/n/2An + yn/2An + C/33/2A33 + V33/2A33-^^33/A33, (30) 
c = V2Wn/Au. (31) 

The eigenvalues are ^(a + 6 it + — 2a6 + 4c^). There are two possible interpretations. 
Either the energies are close together and remain below the isospin 1 states of higher spin, in 
which case we predict two close 1"*" states (which have experimentally not been resolved), or 
their energies are well separated, in which case we predict the observed 1"*" state and a higher 
excited l"*" state that has not yet been seen. From our numerical values for a, b and c we 
calculate the energies to be 0.00194 and 0.00207, which are in fact close together. So we predict 
that the single observed 1'*' isotriplet of states is really an unresolved doublet of isotriplets. 

This matrix diagonalization method is also used to calculate the energy eigenvalues for the 
other values of J'^ and / for which there is more than one allowed state. It is found that in all 
cases, the off-diagonal elements of the matrices analogous to ()28p are small, of the order 10"'^ 
times the diagonal elements. We may therefore consistently assign values of IL3I and \Ks\ to 
each of our calculated energy eigenvalues, as the 'mixing' of states with the same values of J'^ 
and /, but different IL3I and |-ftr3| values, is minimal. The quantum energies of the spin 2 states 
with isospin 1 are: 

^2-,l,2,l = f/ll/An+2;733/A33 + Fil/2Aii+y33/2A33 + 2H^33/A33, (32) 

^2+,i,i,i = 5C/ii/2An + C/33/2A33 + yii/2An + y33/2A33-H^33/A33. (33) 

We calculate that -^2-, 1,2,1 find -^2+, 1,1,1 are 0.00218 and 0.00228 respectively, and so the 2+ 
isotriplet lies above the 2~ isotriplet. Experimentally, however, the 2"*" isotriplet is observed 
below the 2~ isotriplet. 

The quantum energies of the isospin 2 states with spins and 1 are: 

Eo+,2,ofi = 3yii/Aii, (34) 

^l-,2,l,2 = f/ll/2Aii + f/33/2A33 + T4l/Aii+2V33/A33 + 2PF33/A33, (35) 
^1+,2,1,1 = t/ll/2Aii + C/33/2A33 + 5yii/2Aii + y33/2A33-VF33/A33. (36) 

We calculate these values to be 0.00569, 0.00541 and 0.00574 respectively, and so the 1~ state 
lies below the 0^ state, and the 1^ state lies above the 0^ state. Between the O"*" and 1^ states 
there is a further 2"'" state. Experimentally the 0^ isospin 2 quintet includes the ground states 
of beryllium-12 and oxygen-12; low-energy excited 1~ and 2^ states of beryllium-12 are also 
observed. The energy levels of the quantized B = 12 Skyrmion are listed in Table 7. The 
experimental spectrum is in Fig. 9. 



7.1 Three cube interpretation 

In Ref. [16] we estimated the moments of inertia of the B = 8 Skyrmion by treating it as a 
'double cube' of two cubic B = 4 Skyrmions, separated a certain distance along a common C4 
axis, with the cubes rotated around this axis by tt/2 relative to each other. This enabled us 
to estimate the energies of the Skyrmion's quantum states. These estimates agree well with 
our results using the exact B = 8 solution. In this section we apply a similar procedure to the 
B = 12 Skyrmion. 
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=12) 


_ ,.(B=12) 
— '^22 




=12) 


= Mt'' 




=12) 


_ .AB=\2) 
— ^22 


y{B- 
•^33 


=12) 


= 3<=^) 



We work with three B = 4 cubes arranged in an equilateral triangle, meeting at a common 
edge. Each cube is related to its neighbour by a spatial rotation by 2tt/3 followed by an 
isorotation by 27r/3. The isorotation cyclically permutes the values of the pion fields on the 
faces of the cubes, so that these values match on touching faces. We denote by d the distance 
between the centre of each cube and the centre of the triangle. In a similar manner to the case 
of the B = 8 double cube, we use the parallel axis theorem to make estimates for the moments 
of inertia of the -Da/i-symmetric B = 12 Skyrmion in terms of the moments of inertia of the 
B = 4 Skyrmion and the separation parameter d. We use the quadratic interpolation method 
to obtain values for the B = A inertia tensors with m, = 0.685, the optimised value of m for the 
B = 12 sector of the model (see Table 1), to obtain the estimates 

= 3?7jf=^^ = 540, (37) 

= 645 , (38) 

= ?,vi^="^ + ^-M^d\ (39) 

+ 3MAd^ . (40) 

The interpolated value of Ma for m = 0.685 is 589.3. The value of d is chosen using a least 
squares method so that our approximation to Vij is closest to that of the exact B = 12 solution. 
This yields a value for d of 1.92, and hence the estimates 

y/f=i2) = 5756, (41) 
^3(3^=12) = 9014. (42) 

In this simplified picture, Wij vanishes, as it vanishes for the B = A cube. Comparing these 
numbers to those obtained from the exact solution for B = 12 (given in Appendix A. 5), we see 
that the three cube approach has provided good estimates of the inertia tensors Uij and Vij, 
and the inequalities that are satisfied by their elements are right. For the exact solution, Wij 
is non-zero but small. Also, for the exact solution the ratio of Uu to C/33 is closer to 1 than for 
the B = 4 cube. This is because the triangular arrangement of the three cubes is closer to the 
Skyrme crystal, for which Uu = C/33. The accuracy of this approximate inertia tensor shows 
that the Skyrme model is consistent with the intrinsic shape of carbon-12 being an equilateral 
triangle of three alpha-particles. 
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The assumption that Wij vanishes simphfies the quantum Hamiltonian to the sum of a 
symmetric top in space and a symmetric top in isospace: 

IL3I and l-fTsI become good quantum numbers, and the expressions for a selection of the quantum 
energies simphfy to: 
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Their numerical values (using Uu = 540, C/33 = 645, Vu = 5756 and V33 = 9014) are listed in 
Table 7 alongside the corresponding values using the exact solution. We also present the values 
in physical units, obtained using the B = 12 parameter set given in Table 1. 

7.2 Comparison with experimental data 

The ground state of carbon-12 is a O"*" state with isospin 0. It has excited 2"*", 3~ and 4"*" states 
with excitation energies 4.4 MeV, 9.6 MeV and 14.1 MeV respectively. In addition, there may 
be a 4~ state at 13.4 MeV (reassigned from the 2~ state at this energy). As can be seen from 
Table 7, we predict precisely these states, and in the same order. Our predictions for their 
excitation energies are 3.2 MeV, 4.7 MeV, 10.6 MeV and 8.9 MeV respectively. 

Carbon-12 has an excited 0"*" state at 7.7 MeV, the famous Hoyle state. Unfortunately our 
method of rigid-body quantization prohibits two independent spin 0, isospin states. An ex- 
tension of the model might allow this. Perhaps the lowest-lying quantum state of an alternative 
B = 12 solution, such as the C3^-symmetric solution, or a solution with three B = 4 Skyrmions 
in a linear chain, both discussed in Ref. [^, could be interpreted as this excited 0"*" state. 
Additional excited states of carbon-12 are seen experimentally but are not yet predicted in our 
model: for example, a 1^ state at 12.7 MeV. 

We predict two non-degenerate 1"*" states with isospin 1. An isotriplet with J'^ = 1~^ is 
observed, which includes the ground states of boron-12 and nitrogen-12, with average excitation 
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energy 15.1 MeV, to be compared to our value of 12.1 MeV. As mentioned previously, this may 
be an unresolved doublet of isotriplets. However, a number of higher l"*" states with isospin 1 
are seen in the nuclear spectra. We also predict a 2"*" and a 2~ isotriplet. Both of these are 
seen experimentally, but in the opposite order: we predict the 2~ isotriplet to lie below the 
2"*" isotriplet. Experimentally, the 2'^ isotriplet has an average excitation energy of 16.1 MeV, 
compared with 16.5 MeV for the 2~ isotriplet. Higher excited 1~ and 0"^ isotriplets are seen 
experimentally, but we do not predict them in our model. We find two allowed 3"^ states with 
isospin 1, at 16.8 MeV and 18.0 MeV. One such state has been seen in the spectrum of boron- 
12, at 20.8 MeV. An incomplete 3^ isotriplet is observed at 18.5 MeV, which we predict at 
16.4 MeV. We predict a second 3~ state with isospin 1 at 16.8 MeV. States with these quantum 
numbers are seen in the spectra of boron-12, carbon-12 and nitrogen-12, at 20.9 MeV, 20.6 MeV 
and 20.4 MeV, respectively. A further excited 4~ state of boron-12 with isospin 1 exists at 
19.7 MeV, to be compared with our value of 19.0 MeV. 

An (incomplete) J'^ = 0^, isospin 2 quintet, which includes the ground states of beryllium- 12 
and oxygen-12, is observed experimentally with an average excitation energy of 27.7 MeV. We 
predict such an isoquintet at 35.4 MeV. We also predict the existence of a 1~ isoquintet with an 
excitation energy less than that of the 0^ isoquintet. Such a 1~ state is observed in beryllium- 
12, at 2.7 MeV above the O"*" ground state. An (incomplete) 2"*" isoquintet is experimentally 
observed with an average excitation energy roughly 2 MeV above the 0+ isoquintet. In our 
model three 2"*" states with isospin 2 are allowed, just above the 0^ state. 

In summary: The model describes the energy spectra of nuclei of mass number 12 especially 
well. The rotational band of carbon-12 is very well reproduced, along with some of the exper- 
imentally observed isospin 1 triplets and isospin 2 quintets. However the observed isotriplets 
with J'" = 0"*" and 1~ do not appear as quantum states of our L'3/i-symmetric Skyrmion. Nei- 
ther does the Hoyle state. They may arise from the quantization of further modes, or appear as 
quantum states of an alternative Skyrmion. The molecular rotational band of beryllium-12 in 
the range 10-20 MeV above the beryllium-12 ground state [27] may also be explained in terms 
of an alternative Skyrmion. 
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Figure 9: Energy level diagram for nuclei of mass number 12. Mass splittings between nuclei 
are adjusted to eliminate the proton/neutron mass difference and remove Coulomb effects, as 
described in Ref. l23l . 
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Wave function ^ 
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Table 6: FR-allowed quantum states of the B = 12 Skyrmion (normalised). 
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Table 7: Energy levels of the _B = 12 Skyrmion, using both the exact solution and the three cube 
interpretation. To each of the quantum states of the exact solution there correspond dominant 
values of IL3I and \K^\. 
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8 Conclusion 



The Skyrme model provides a model of nuclei which unifies spin and isospin as collective excita- 
tions of topological solitons. We have been able to reproduce the energy spectra of a number of 
light nuclei to a good degree of accuracy, by quantizing Skyrmions as rigid bodies in space and 
isospace, and by parametrising the model separately for each mass number. The model is fitted 
to nuclear charge radii and masses. Remarkably, this leads to good predictions for moments of 
inertia, spin splittings and isobar splittings. The quantum states we have calculated agree well 
with what is experimentally observed, with the correct spins and parities. Among the successes 
of the model is the prediction of the existence and the excitation energies of the rotational 
bands of beryllium-8 and carbon-12. We predict that neither beryllium-10 nor carbon-10 has a 
l"*" state, in agreement with experiment. Our predictions of the quantum states of the nuclear 
isobars with non-zero isospin is quantitatively good, and the spin splittings are quite good. 
However, the spin states do not always appear in the correct order. A few, so far unobserved 
states have been predicted, notably some negative parity states of lithium-8 and boron-8, and 
some high spin states including a 4"*" state of helium-4 and a 6^ state of carbon-12. Several 
observed states are not predicted by the rigid-body quantization of Skyrmions. To understand 
these one will need to consider further Skyrmion degrees of freedom. 

One can be pleased with the general trend of isospin excitations. In each case, the / = 0, 1 
splitting is of the order of 10 MeV, the / = 0, 2 splitting of the order of 20-30 MeV, and the 
/ = 0, 3 splitting of the order of 60 MeV. Isospin splittings for spherically averaged Skyrmions, 
over a wider range of baryon numbers and isospins, have been estimated in Ref. [13j. The 
Skyrme model predicts that the isobar splittings rise in proportion to /(/ -|- 1). Curiously the 
experimental data do not rise quite as fast, but our prediction for the coefficient of /(/ -|- 1) is 
roughly correct. 

It would have been possible to have chosen just one parameter set, as in previous work [16] . 
We could have chosen m close to 1 , the energy scale close to 6 MeV, and the length scale close 
to 1.6 fm, for example. This would have provided good qualitative results but would not have 
given such a good fit to charge radii and energy splittings. 

Further work is needed on the electromagnetic form factors and transition amplitudes of 
light nuclei within the Skyrme model, as these provide more information about the internal 
structure of nuclei. It would also be desirable to consider the effect of vibrational modes and 
the break-up of Skyrmions into clusters, for example, modelling the break-up of lithium-6 into 
helium-4 and the deuteron. Finally, we would like to extend this analysis to Skyrmions and 
nuclei with B = 7, 9, 11 and beyond 12. 
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A Numerical Results 



The static energies, mean charge radii and the non-zero elements of the inertia tensors for 
B = A, 6, 8, 10 and 12 have been calculated for m = 0.5, 1 and 1.5, and are shown below. The 
interpolated vahics for the preferred values of m are shown in the final columns. An estimate of 
the numerical errors can be made by calculating the off-diagonal inertia tensor elements, which 
should be identically zero in each case. The calculated ratio of off-diagonal elements to non-zero 
diagonal elements is of the order of 10"^ or less. 
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A.4 5 = 10 
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A. 6 Quadratic interpolation between three points 

To obtain approximations to the inertia tensors, static energies and mean charge radii for a 
given value of m, we use the method of quadratic interpolation between three points. We know 
the values of a property p at m = 0.5, 1 and 1.5, and we make the ansatz 

p{m) = «! + a2m + a3vr? . (45) 

Let p = (p(0.5),p(l),p(1.5))"^. The vector a. = (ai,Q;2,a3)^ is obtained by inverting the 
expression 

M • ct = p , (46) 

where 

1 0.5 0.25 

M = ( 1 1 1 I , and therefore M"^ = | -5 8 -3 | . (47) 
1 1.5 2.25 

The interpolated values are given in the right-hand columns of the tables in Appendices A.l to 
A.5. 
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